Entanglement signature in the mode structure of a single photon 
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It is shown that entanglement, which is a quantum correlation property of at least two subsystems, 
is imprinted in the mode structure of a single photon. The photon, which is emitted by two coupled 
cavities, carries the information on the concurrence of the two intracavity fields. This can be useful 
for recording the entanglement dynamics of two cavity fields and for entanglement transfer. 
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An atom interacting with a quantized radiation-field 
mode in a high-Q optical cavity plays an important role 
in quantum optics, for a review see, e.g., Ref. [H. The 
ability to coherently control individual quantum system, 
and in particular the quantum control of single-photon 
emission from an atom in a cavity, is a key requirement 
in various applications of quantum networks for distribu- 
tion and processing of quantum information [2, Ji, 4, S,^. 
Recently, single-photon sources operating on the basis of 
adiabatic passage with just one atom trapped in a high-Q 
optical cavity have been realized 0, H, [^.llOj] . In this way, 
the adjustment of the spatiotemporal profile of single- 
photon pulses has been achieved [ll|, [ij] ■ Moreover, the 
generation of single photons of known circular polariza- 
tion emitted into a well-defined spatiotemporal mode has 
been possible [13j , and an atom-photon quantum inter- 
face involving atom-photon entanglement has been real- 
ized fl^. More recently, the amplitude modulation in 
the photon emission on a single atom-cavity system has 
been studied theoretically [l5| and experimentally [iGj . 
In addition, photon-photon entanglement with a single 
trapped atom in a high-finesse optical cavity has been 
performed [l7j . 

In the present contribution, in view of the widespread 
applications of cavity-assisted single-photon sources, we 
study single-photon emission from a system consisting of 
two coupled atom-cavity subsystems in a cascaded con- 
figuration [3, [3 ■ The mode structure of the radiated 
photon strongly depends on the entanglement between 
the two intracavity fields and it sensitively depends on 
the presence or absence of an atom in the second cavity. 
We show how the entanglement of the intracavity fields 
can be experimentally determined. 

The system under study consists of two atom-cavity 
subsystems A and B, where the source subsystem A is 
cascaded with the target subsystem B, cf. Fig. [T] The 
cavities have three perfectly refiecting mirrors and one 
mirror with transmission coefficient T <C 1. In the two 
subsystems A and B we consider a two-level atomic tran- 
sition of frequency ujk (related to the atomic energy eigen- 
states |lfc) and |Ofc)) coupled to a cavity mode of fre- 
quency io'f., where k ~ a,b denotes the subsystem. The 
cavity mode is detuned by Afc from the two-level atomic 
transition frequency, ujk = oj'j, + A^, and is damped by 
losses through the partially transmitting cavity mirrors. 
In addition to the wanted outcoupling of the field, the 
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FIG. 1: The cascaded system consists of two atom-cavity sub- 
systems A and B. A photodetector PD monitors the radiation 
field. 



photon can be spontaneously emitted out the side of the 
cavity into modes other than the one which is preferen- 
tially coupled to the resonator. Moreover, the photon 
may be absorbed or scattered by the cavity mirrors. 

To describe the dynamics of the system we use the fol- 
lowing master equation for the reduced density operator 
p{t) of the system: 



dp{t) 
dt 



1 

ih 



The Hamiltonian is given by 



H = Ha + Hb+ ih 



(1) 



(2) 



where Ha and Hb describe the atom-cavity interaction 
in the two subsystems A and B, respectively, and, in the 
rotating-wave approximation, are given by 

Ha = hga (aiio + aUoi) + /lAain , (3) 



and 



Hb - %6 ( SSio + StSoi ) + h^bBll . (4) 



2 



The third term in E g. pt describes the coupHng between 
the two cavities [Tsl, Il9l|. In these expressions, a (a^) and 
b (6^) are the annihilation (creation) operators for the 
cavity fields A and B, respectively. We have also defined 
Aj = \ia){ja\ {hj = 0, 1), and B,j = \ib){jb\ («, J = 0, 1). 
In addition, gk is the atom-cavity coupling constant and 
Kfc the cavity bandwidth, and the phase (f> is related to 
the spatial separation between the source and the target, 
cf. [20|. The jump operators Ji are defined by 

ii = y/i^a + y/K^e^'-'^b , (5) 

which describes photon emission by the cavities; 



J2 



(6) 



are associated with photon absorption or scattering by 
the cavity mirrors; and 



ra^oi , 



ThBa 



(7) 



are related to spontaneous emission by the atoms. Here 
kJ. and Ffc are the cavity mirrors' absorption (or scat- 
tering) rate and the spontaneous emission rate of the 
two-level atom, respectively. Note that the operator Ji 
contains the superposition of the two fields radiated by 
the two cavities, due to the fact that the radiated photon 
cannot be associated with photon emission from either A 
or B separately. 

To evaluate the time evolution of the system we use a 
quantum trajectory approach [20l.[2ll.[2^. Let us consider 
the system prepared at time to = in the state |a) = 
|1, 0,0,0), which denotes the atom A in the state |la), 
the cavity A in the vacuum state, the atom B in the state 
I Ob) , and the cavity B in the vacuum state. Similarly, we 
define \b) = |0, 1,0,0), |c) = |0, 0,1,0), \d) = |0, 0,0,1), 
and |e) = |0, 0,0,0). To determine the state vector of 
the system at a later time i, assuming that no jump has 
occurred between time to and t, we have to solve the 
nonunitary Schrodinger equation 



ih-\,Pnoit)) = H' l^b^^it)) , (8) 

where H' is the non-Hermitian Hamiltonian given by 



H' ^ H 



where we have defined Ka — Ka + n'a and — Kh + K^,. If 
no jump has occurred between time to and t, the system 
evolves via Eq. ([8]) into the unnormalized state 

l^no(t)) = a{t)\a) -f m\b) + l{t)\c) + 5{t)\d) . (10) 

The evolution governed by the nonunitary Schrodinger 
equation ([5]) is randomly interrupted by one of the five 



1=1 
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kinds of jumps Ji, cf. Eqs. ©-([7]). If a jump has occurred 
at time tj, tj G ito,t\, the wave vector is found collapsed 
into the state |e) due to the action of one of the jump 
operators 



id^„o(tj)) |e) (i = l,...,5). 



(11) 



In the problem under study we may have only one jump. 
Once the system collapses into the state |e), the nonuni- 
tary Schrodinger equation ([8]) lets it remain unchanged. 
The density operator p(t) is then obtained by perform- 
ing an ensemble average over the different trajectories at 
time t, yielding the statistical mixture 



m = |^„o(i))(Vi„o(t)| + k(t)ne)(e| 



(12) 



where |e(t)p = 1 - (V'„o(t)|</'„o(i)>- The values |a(t)p, 
l7(^)P> and |e(i)|2 represent the probabil- 

ities that at time t the system can be found either in |a), 
1 6), |c), \d), or |e), respectively. 

In order to determine a{t), I3{i), j{t), and S{t), we 
have to solve the nonunitary Schrodinger equation, cf. 
Eqs. (HI) and ([9]), which leads to the inhomogeneous sys- 
tem of differential equations 

a{t) = ~i (A, - iTa/2) a{t) - igaP{t) , 
/3(t) = -igaa{t) ~ {Kj2)(3{t) , 
7(t) = -z (Ab - tTb/2) 7(t) - igbS{t) , 
5{t) = -igblit) - {Kbl2)5{t) - ^WJTbe'* I3{t) . 

(13) 

For the initial conditions a(0) = l, /3(0) = 0, 7(0) = 0, and 
(5(0) = 0, and defining 



^^-Agl-iKkUk-i^- 



- A.-? 



we get, similarly as done in [23|, the solutions 

'Ka/2-i{^a- iTa/2) 



(14) 



a{t) = 

m = 



■ sinh 



^-[(K^+Ta.)/i+iA^/2]t 



cosh(^^ 

Sinhf — ) f,-[(Ka+T^)/i+^AJ2]t ^ 
\ ^ / 

9b {/+ it) [5_ {t) + h+ {t)] - /_ it) [g+ {t) + h^ m 



Kb-Tb Ab 

1 

4 2 

Kb-Tb .Afc 
4 * 2 ^ 



2 
2 



Ut)[g+{t) + h^it)] 
f+it)[g^it) + h+{t)]. (15) 



Here we have defined, 

^^^^^ ^ gaV^a^bg"^ ^f-(g,+r,)/4-»At,/2±a,/2]t 



9±{t) 



g[(0„±f2b)/2-T-zA]t _ I 

{na±^b)/2-T -iA ' 



(16) 



(17) 



3 



and 



h±{t) 



g-[(f2„±06)/2+T+iA]t _ I 
(f7a±f^fc)/2 + T + iA 



(18) 



where T^{Ka-Kb+Ta-Tb)/4: and A=(A,-Ab)/2. In 
the case of equal parameters for the two subsystems A 
and B, the solutions for j{t) and 6{t) simplify to 
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tsinh 



-[(K+r)/4+iA/2]t 
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(19) 



where k — Kfc, K = Ka = -R'h, 
Tfe, g^ga^gt, and f7 = f2a=:fi{,. 

In the system under study, because only one atom 
is initially excited, the two intracavity fields constitute 
a pair of entangled qubits, for a detailed discussion of 
single-particle entanglement, see [13]. An appropriate 
measure of the entanglement for a two-qubit system is 
the concurrence [25| . To derive an expression for the con- 
currence between the two intracavity fields we consider 
the density operator obtained by tracing over the atomic 
states for the two subsystems . Pr ,m(t) = Trat [/5(i)]. It is 
easy to show, following Ref. [2^, that the concurrence 
between the two intracavity fields is given by 



c[p,,,{t)] = 2\pit)\m\ 



(20) 



Note that for equal parameters for the two subsys- 
tems, and for g 3> i^, F, A, the concurrence is given by 
C[pcav(i)]=^«isin!(fft)e[-(^+n*/2]. 

Following [2^, [23 , we consider a photon in the mode 
^i, the mode escaping from the cavities and going to the 
photodiode PD. It is described by the normalized func- 
tion ^i{t) of amplitude envelope and phase (f>i(t), 

= with 

dtMt)?^ / dtC^{t)^l. (21) 
/o Jo 

When a photon is in the mode ^j, whose amplitude en- 
velope C,i{t) does not change significantly in the detec- 
tion time resolution T, the response probability of the 
detector of quantum efficiency ri within a time interval 
[t - T/2, t + r/2] is given by [l| 

/'D(t) =?7Prad(00)Cf(i)r. (22) 

Here Prad(oo) — linit_+ooPrad(Oi where the function 
Pradit) represents the probability that a photon is ra- 
diated by the cascaded system in the time interval [0, i], 
which reads as 



Prad(i)=/ dt'{JlJi)t'^Kal dt' \(3{t')\^ + K^J ^ dt' \S{t')\' 



+ 2y/J^ / dt'Re[f3*{t')S{t')t 



(23) 
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FIG. 2: The function 2Re[/3* (t)5(t)e-''^] is shown for equal 
parameters for the two subsystems A and B, for g/K — 5, 
k/K = 0.9, A/A' = 0.1, r/K = 0.2 (solid line), when Eq. ^ 
applies. The case when no atom is present in the second 
cavity, i.e. gb — 0, is also shown (dashed line). 



Since S{t) contains an overall factor e"^, cf. Eqs. (fT5|) and 
the phase (j) is irrelevant in Eq. 
The probability to measure between time t — T/2 and 
t+T/2 a "click" at the detector is equal to the probability 
to have a jump Ji in the same time interval, so that using 
Eq. (O, we get 



(24) 



Puit) = 77Tr [mjlJi] T = vT[Ka\m\ 
+ 2^]^Re[l3*{t)S{t)e~"^]Y 
Comparing this with Eq. ([22]) we obtain 
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Ctit) 



Prad(oo) 



(25) 

Note that Eq. ^ is correctly fulfilled. 

Let us now analyze in more details the term 
2Re[/3*(i)5(i)e-''^] in Eq. Writing P{t) = 

|/3(i)|e*"^'^(*' and S{t) = e''^|J(i)|e*"^^(*), yields 

2Re[l3*it)S{t)e-''^]^C[p,,^it)]cos[Mt) - Mi)] > (26) 

where C[pcnvit)] is the concurrence between the two in- 
tracavity fields, cf. Eq. ([20| . In this respect, Eq. ([25]) 
clearly shows that the mode structure of the radiated field 
depends not only on the two intracavity fields, i.e. |/3(i)P 
and |iy(t)P, but also on the entanglement established be- 
tween them. This represents an interference between the 
possibility to have the photon in one or in the other cav- 
ity. For equal parameters for the two subsystems, and 
for g ^ K, F, A, one obtains the relation 



2Re[(3*{mt)e-'^] ~ -C[p,,^{t)] 



(27) 



In this case the concurrence can be experimentally de- 
rived by using the combination of two measurements. 
The first one, by using only cavity A gives |/3(t)|, via the 
relation P^(t) = r]KT\P{t)\'^ , cf. The second mea- 

surement, by using both cavities, gives \5{t)\ via the re- 
lation PY,{t)/P^{t) = (1- |5(0|/|/3(i)l)'- Knowing \(3{t)\ 
and \5{t)\, the concurrence is obtained from Eq. (|^D|) . 
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FIG. 3: The amplitude envelope C»(*)\/Prad(oo)/K for the 
mode of the radiated field is shown for equal parameters for 
the two subsystems A and B, for g/K = 5, k/K — 0.9, 
A/K = 0.1, r/K = 0.2 (full line), and for no atom in the 
second cavity (dashed line) . The case where the subsystem B 
is absent, i.e. for Kt — 0, is also shown (dotted line). 

In Fig. m we show the term 2Re[P*{t)dit)e-"l'] under 
conditions when it represents the negative concurrence 
according to Eq. ((27)) . We also show the case when no 
atom is present in the second cavity. In both cases the 
entanglement between the two intracavity fields gives a 
significant contribution to the mode structure of the ra- 
diated photon, cf. Eq. ([25|l . 



The amplitude envelope for the mode of the radiated 
field is shown in Fig. [31 for equal parameters of the two 
subsystems and for the case when no atom is present in 
the second cavity. The case when the subsystem B is ab- 
sent, i.e. for Kb — 0, is also shown, reproducing the result 
obtained in [15]. The shown mode structures earring the 
entanglement signature could be realized and observed 
by extending the experimental setup described in p^ . 
By measuring the arrival time distribution of the photon 
radiated from a system with equal cavity parameters, one 
may determine the full dynamics of the concurrence and 
hence the entanglement dynamics of the two intracav- 
ity fields in the strong coupling regime. This regime has 
been realized in recent experiments [3, ■ 

In conclusion, the dynamics of a system consisting of 
two atom-cavity subsystems has been analyzed under re- 
alistic conditions with losses. For properly chosen param- 
eters, the mode function of the single photon escaping 
from the cavities reflects the full dynamics of the concur- 
rence of the two intracavity fields, while they continue to 
interact with the two atoms. This allows one to detect 
the entanglement dynamics of two cavity fields, and may 
be useful for transferring the information on entangle- 
ment by a single photon over a large distance. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft. 



[1] S. Haroche and J.-M. Raimond, Exploring the Quantum 

(Oxford University Press, Oxford, 2006). 
[2] J.I. Cirac, P. ZoUer, H.J. Kimble, and H. Mabuchi, Phys. 

Rev. Lett. 78, 3221 (1997). 
[3] S. Brattke, B.T.H. Varcoe, and H. Wahher, Phys. Rev. 

Lett. 86, 3534 (2001). 
[4] E. Knill, R. Laflamme, and G.J. Milburn, Nature 409, 

46 (2001). 

[5] C. Monroe, Nature 416, 238 (2002). 

[6] H.J. Kimble, Nature (London) 453, 1023 (2008). 

[7] A.S. Parkins, P. Marte, P. ZoUer, and H.J. Kimble, Phys. 

Rev. Lett. 71, 3095 (1993). 
[8] M. Hennrich, T. Legero, A. Kuhn, and G. Rempe, Phys. 

Rev. Lett. 85, 4872 (2000). 
[9] J. McKeever, A. Boca, A.D. Boozer, R. Miller, J.R. Buck, 
A. Kuzmich, and H.J. Kimble, Science 303, 1992 (2004). 
[10] M. Hijlkema, B. Weber, H.P. Specht, S.C. Webster, 
A. Kuhn, and G. Rempe, Nature Physics 3, 253 (2007). 
[11] A. Kuhn, M. Hennrich, and G. Rempe, Phys. Rev. Lett. 

89, 067901 (2002). 
[12] M. Keller, B. Lange, K. Hayasaka, W. Lange, and 

H. Wahher, Nature (London) 431, 1075 (2004). 
[13] T. Wilk, S.C. Webster, H.P. Specht, G. Rempe, and 

A. Kuhn, Phys. Rev. Lett. 98, 063601 (2007). 
[14] T. Wilk, S.C. Webster, A. Kuhn, and G. Rempe, Science 

317, 488 (2007). 
[15] C. Di Fidio, W. Vogel, M. Khanbekyan, and D.- 



G. Welsch, Phys. Rev. A 77, 043822 (2008). 

[16] J. Bochmann, M. Miicke, G. Langfahl-Klabes, C. Erbel, 

B. Weber, H.P. Specht, D.L. Moehring, and G. Rempe, 

Phys. Rev. Lett. 101, 223601 (2008). 
[17] B. Weber, H.P. Specht, T. Miiller, J. Bochmann, 

M. Miicke, D.L. Moehring, and G. Rempe, Phys. Rev. 

Lett. 102, 030501 (2009). 
[18] H.J. Carmichael, Phys. Rev. Lett. 70, 2273 (1993). 
[19] C.W. Gardiner, Phys. Rev. Lett. 70, 2269 (1993). 
[20] H.J. Carmichael, Statistical Methods in Quantum Optics 

2, (Springer, Berhn, 2008). 
[21] J. Dahbard, Y. Castin, and K. M0lmer, Phys. Rev. Lett. 

68, 580 (1992). 

[22] R. Dum, A.S. Parkins, P. ZoUer, and C.W. Gardiner, 

Phys. Rev. A 46, 4382 (1992). 
[23] C. Di Fidio and W. Vogel, Phys. Rev. A 78, 032334 

(2008). 

[24] S.J. van Enk, Phys. Rev. A 72, 064306 (2005). 

[25] W.K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[26] K.J. Blow, R. Loudon, S.J.D. Phoenix, and T.J. Shep- 
herd, Phys. Rev. A 42, 4102 (1990). 

[27] T. Legero, T. Wilk, A. Kuhn, and G. Rempe, Appl. Phys. 
B 77, 797 (2003); Adv. At. Mol. Opt. Phys. 53, 253 
(2006). 

[28] C.J. Hood, T.W. Lynn, A.C. Doherty, A.S. Parkins, and 

H. J. Kimble, Science 287, 1447 (2000). 



